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We study how well the Gaussian approximation is valid for computing the covariance matrices of 
the convergence power and bispectrum in weak gravitational lensing analyses. We focus on its impact 
on the cosmological parameter estimations by comparing the results with and without non-Gaussian 
error contribution in the covariance matrix. We numerically derive the covariance matrix as well 
as the cosmology dependence of the spectra from a large set of A^'-body simulations performed for 
various cosmologies and carry out Fisher matrix forecasts for tomographic weak lensing surveys with 
three source redshifts. After showing the consistency of the power and bispectra measured from our 
simulations with the state-of-the-art fitting formulae, we investigate the covariance matrix assuming 
a typical on-going survey across 1,500 deg^ with the mean source number density of 30 arcmin"^ 
at the mean redshift Zs = 1.0. Although the shape noise contributes a significant fraction to the 
total error budget and it mitigates the impact of the non-Gaussian error for this source number 
density, we find that the non-Gaussian error degrades the cumulative signal-to-noise ratio up to the 
maximum multipole of 2,000 by a factor of about 2 (3) in the power (bi-) spectrum analysis. Its 
impact on the final cosmological parameter forecast with 6 parameters can be as large as 15% in the 
size of the one-dimensional statistical error. This can be a problem in future wide and deep weak 
lensing surveys for precision cosmology. We also show how much the dark energy figure-of-merit is 
affected by the non-Gaussian error contribution and demonstrate an optimal survey design with a 
fixed observational time. 

PACS numbers: 98.80.Es 



I. INTRODUCTION 



Cosmological weak gravitational lensing has been be- 
coming the focus of attention as a poviferful tool to probe 
the distribution of the matter in the universe, since its 
first detections [ll-lll • The lights from distant galaxies are 
bent by the gravitational potential of intervening large 
scale structures, which generate coherent deformation of 
galaxy images and this is the so-called cosmic shear. We 
can directly see the distribution of matter, and measure 
the statistical quantities such as the power spectrum and 
bispectrum of mass fiuctuations on cosmological scales 
by analyzing this coherent pattern. 

Weak lensing can also be a powerful probe of the nature 
of dark energy. The growth rate of mass clustering can 
be measured by lensing tomography [e.g.,[5l-Q which in 
turn provides tight constraints on the equation of state 
of dark energy. For this purpose, a number of ambi- 
tious wide-field surveys have been proposed, such as Sub- 
aru Hyper Suprime-Cam Survey (HSC^) [1], the Dark 
Energy Survey (DES^) @, the Large Synoptic Survey 
Telescope (LSST^) the Wide-Field Infrared Survey 



Telescope (WFIRST^), and Euchd^ [TTl. 

Unlike the cosmic microwave background, the distribu- 
tion of matter in the present universe which determines 
the weak-lensing convergence field is highly non-linear 
and non-Gaussian, reflecting the non-line ar p rocesses 
that accompanied structure formation |la - [l4j . Thus, 
some of the cosmological information initially stored in 
the two-point correlation function and/or the power spec- 
trum when the density field was in the linear stage is no 
longer present in two-point statistics of the non-linear 
convergence field. If we want to draw as much infor- 
mation as possible from the convergence field, it is of 
great importance to add the information contained in 
the higher order statistics such as the bispectrum on top 
of the two-point statistics or to resort to alternative non- 
standard methods. 

Recently, it was suggested that the two-point statistics 
of the logarithmic transformed non-linear weak lensing 
field may contain more information than the two-point 
statistics of the original field before transformation [e.g., 
[l3| . Seo et al. [l^ showed that, by Taylor-expanding 
the log-transformed field, most of the improvement by 
the transformation can be explained by the information 
originally contained in the bispectrum, suggesting that 
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the log-transform is a way to draw some information in 
the higher-order statistics back to the two-point statis- 
tics. However, it was also shown that the log-transform 
is advantageous when one is interested in a single param- 
eter such as the amplitude of the power spectrum. Seo 
et al. [Til found that there is only little improvement in 
the constraints on multiple cosmological parameters af- 
ter log-transform in the presence of shape noise that are 
expected from future weak lensing surveys. 

We therefore focus on the weak lensing bispectrum it- 
self as a higher-order statistic in this study. We study 
the usefulness and complementarity of the lensing bispec- 
trum compared to the power spectrum using 1,000 ray- 
tracing weak-lensing maps generated in Seo et al. [15|. 
We consider all the possible triangular configurations 
available from a given range of multipoles and examine 
the impact of the bispectrum on cosmological parame- 
ter estimation taking non-Gaussian covariance matrices 
obtained from ray-tracing simulations into account. In 
doing this, we perform a Fisher matrix analysis of a tomo- 
graphic survey with three different source redshifts using 
the fully non-Gaussian covariance matrix and derivatives 
of the spectra wrt the cosmological parameters directly 
measured from the simulations. 

The paper of Kayo et al. [3 recently appeared is based 
on a similar analysis. They mainly focused on the co- 
variance matrix itself and examined the impact of the 
new contribution coming from the number fluctuation 
of massive halos in a finite survey area based on halo 
model [l3l, which we hereafter refer to as halo sample 
variance (HSV). In this paper, we will extend their anal- 
ysis to the information content of the power and the bis- 
pectrum tomography in terms of the cosmological param- 
eter constraints. 

This paper is organized as follows. In Section |lTl we 
briefly review the basic theoretical expectations of the 
lensing power and bispectrum and their covariances. We 
describe the details of our A^-body simulations used in 
this paper and data analysis in Section lllll In Section HVl 
we present the detailed comparisons between the sim- 
ulation results and previous models for the weak lens- 
ing power and bispectrum. After studying the bispec- 
trum covariance using the simulations in Section |Vl we 
study cumulative signal-to-noise ratio of the power spec- 
tra, bispectra and the joint measurement of them in Sec- 
tion IVll In Section IVlll we present how significant the 
non-Gaussian errors of the weak lensing power and bis- 
pectrum arc in constraining the cosmological parameters 
by using the Fisher matrix analysis. In Section IVllll we 
examine the dark energy Figure-of-Merit in the presence 
of the non-Gaussian corrections to the covariance matrix, 
and demonstrate an optimal survey design by changing 
the mean number density of source galaxies and the sur- 
vey area under the condition that the total observation 
time is fixed. Finally, Section HXl is devoted to conclusion 
and discussion. 



II. PRELIMINARIES 

A. Lensing power spectrum and bispectrum 

In this paper, we consider an ideal source galaxy distri- 
bution where all of them are located at a single redshift 
when we compare with simulations. By using the Born 
approximation, the weak-lensing convergence field can be 
written as a weighted projection of the three-dimensional 
density contrast [e.g.. Iisl. [Tgj 
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where 6 is the two-dimensional vector denoting the angu- 
lar position on the sky, x is the comoving distance, Xs is 
the comoving distance to the source, Jk (x) is the comov- 
ing angular diameter distance, and W{x) is the lensing 
weight function defined as 

mx) = ^^^^^^%^^^^^^(l + .). (2) 

To compute the power and bispectrum of the conver- 
gence, we employ the flat-sky approximation which is 
sufficiently accurate over angular scales of our inter- 
est. Within this approximation, the lensing convergence 
field is decomposed into angular modes based on two- 
dimensional Fourier transform: 



(3) 



For a two-dimensional homogeneous and isotropic ran- 
dom field, one can define the convergence power spectrum 
and bispectrum as 

{k{h)k{l2)) = {27t)^Sd{Ii + l2)P.{li), (4) 

and 

{k{h)~K{h)~Kih)) = {27:)^dD{h+l2+h)B^ill,l2j3), (5) 

where S]j{l) is the Dirac delta function. By using the 
Limber approximation [20l . [2]| . the conver geii ce p ower 
spectrum and bispectrum are given by [e.g.. [iSl [l9j 
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where fc^ = h/ fxix)^ Ps{k]z), and Bs{ki,k2,kz; z) are 
the three-dimensional power spectrum and bispectrum 
of the matter density contrast at redshift z. The non- 
linear gravitational evolution of Ps and Bs significantly 
enhances the amplitude of the lensing power and bispec- 
trum respectively at I ^ 100 for source redshift = 1.0 
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(see, Figures [T] and [2]) . Therefore, we need to take non- 
linear evolution effect into account for weak lensing stud- 
ies. We employ some fitting formulae for the three dimen- 
sional spectra in analytically evaluating the convergence 
spectra; halof it proposed by Smith et al. and a re- 
fined version of that proposed by Takahashi et al. [l^l 
(hereafter revised halof it) for Pg and Scoccimarro 
and Couchman [23| (hereafter SCOl) and Gil-Marfn et al. 
di] (hereafter Gil-Marinl2) for Bs in Sec|lVl The two 
fitting functions for the bispectrum explicitly include Pg, 
and we will use halof it and revised halof it for that. 

The power spectrum measures only the Gaussian infor- 
mation, and one can complement it by relying on some di- 
rect measures of non-Gaussianity such as the bispectrum. 
The practical advantage of the bispectrum is that the in- 
trinsic ellipticities do not contaminate the measurement, 
as long as its intrinsic distribution is symmetric with zero 
mean. The measured bispectrum is therefore an unbiased 
estimator under this condition. The B-modes generated 
from the intrinsic ellipticity alignment or higher order 
lensing effects also do not contaminate the bispectrum 
measurement thanks to the parity invariance [26|. Thus, 
systematic errors associated with bispectrum measure- 
ment are smaller than the power spectrum measurement. 



B. Covariance matrices of the lensing power 
spectrum and bispectrum 



small scales. This additional variance, HSV, is expressed 
as 

CovHsv[P«(0,P.(n]=f^dx(g^)^ 
X / dM^b{M)\kM{l)\' I dM'^b{M')\kM'{lT 
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where d^V^/dxdfi is the comoving volume per unit radial 
comoving distance and unit solid angle, and is given by 
for a flat universe, dn/dM is the ensemble-averaged 
halo mass function, b{M) is the halo bias parameter, and 
k.m{1) is the angular Fourier-transform of the convergence 
field generated by the density profile of a halo with mass 
M . Also, W{x) is the Fourier transform of the survey 
window function. It should be noted that HSV contri- 
bution does not necessarily scale with l/f2s unlike other 
covariance terms because the sample variance depends 
on fig via the shape of the linear power spectrum Pj"(fc). 
The contribution of Equation (|lip arises for any finite- 
volume survey because the halo distribution has modu- 
lations due to the biased density fluctuations over the 
survey window. 

Meanwhile, the covariance matrix of the convergence 
bispectrum is defined as a sum of five terms (igI l28l| : 



The covariance matrix of the convergence power spec- 
trum between Pk{1) and Pk.{1') can be expressed as a sum 
of the Gaussian and non-Gaussian contributions I17ll27l: 



Cov[P«(/),P«(O] = -r^P.(0''5/S 



(8) 

where 6l^[, is the Kronccker delta function, ilg is the sur- 
vey area, and r„ is the lensing trispectrum defined as 

{K{li)fi{l2)K{l3)fi{h)) = {2TTfSD{h234:)T^illj2,l3j4), 

(9) 

where we have introduced a shorthand notation Z1234 = 
h + h + h + h- In the above, Ni denotes the number 
of modes around a bin labeled by I and is approximately 
given by Ni = Asr2s/(27r)2 with = 2ttIAI being the 
area of the two-dimensional shell around that bin. There- 
fore, I denotes the mean radius of the annulus. In the 
Limber approximation, is simply a projection of the 
three-dimensional mass trispectrum Tg given by 
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However, there is an additional contribution to the 
non-Gaussian covariance, which becomes significant on 
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The first term is proportional to the triplet of the lensing 
power spectrum with an amplitude given by the geomet- 
rical factor: 



(27r)^AiJ2,i3,i'i,i2^i3 
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where 
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and 
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A-'ih,l2, h) = -^^2lili + 2lili + 2lili -It^lj- 

(15) 

if 1^1 — /2I < ^3 < h + h is satisfied, else A ^(^1,^2, '3) = 0. 
This factor shows the area of a triangle with side lengths 
h, h, and ^3. We refer to the first term of Equation (fT2|) 
as the Gaussian contribution, while the other four terms 
denote the non-Gaussian contributions that arise from 
the connected three-, four-, six-point function of the con- 
vergence field, and the number fiuctuations of massive 
halos in a finite survey area (see, [l6j for the exact expres- 
sions, while 12 8i missed the final term which is dominant 
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over the other covariance term at I ^ 1000 even if survey 
area is a few thousand (see, Figure 2 in Kayo et al. (l6|). 

Only the first term in Equation ([T2|) is usuahy dis- 
cussed in previous statistical analyses of the cosmolog- 
ical fields in the literature 0, except for [l^, just 
because of simplicity and/or difficulty of calculation of 
the last four terms. We will carefully examine how well 
the approximation of Gaussianity is valid for computing 
the bispectrum covariance and its impact on cosmological 
parameter estimations, by comparing with the fully non- 
linear covariance matrix measured from a large ensemble 
of ray-tracing simulations. 



III. NUMERICAL SIMULATION 

A. Simulation design 

In order to study the impact of the non-Gaussian er- 
ror of the convergence power and bispectrum on cosmo- 
logical parameter estimations, we perform a large set 
of ray-tracing simulations through large-volume, high- 
resolu tion TV-body simulations of structure formation |17l . 

We use a modified version of the Gadget-2 
33[ for the A^-body simulations. The matter den- 
sity fields in quasi-light cone volumes are constructed 
by combining 2 x 200 realizations of iV-body simula- 
tions performed in cubes with 240 and 480/i~^Mpc on 
a side, and we perform ray-tracing simulations through 
these volumes. Wc employ 256'^ particles for each TV- 
body simulation. For our fiducial cosmology, we adopt 
the standard ACDM model with density parameter of 
matter = 0.238, baryon fib = 0.0416, and dark en- 
ergy r^A = 0.762 with the current value of the equa- 
tion of state parameter wq = —1 and its time evolution 
Wq = 0, the primordial spectrum with the spectral in- 
dex Us — 0.958 and the normalization As = 2.35 x 10~^, 
and the Hubble parameter h = 0.732, which are consis- 
tent with the WMAP 3-year results [sil- The amplitude 
of the linear density fluctuations in a sphere of radius 
8/i~"'^Mpc at present time is = 0.759 in this cosmol- 
ogy. We assume three delta-function-likc source redshifts 
at Zg = 0.6, 1.0, and 1.5 to perform a tomographic study. 
Using ray-tracing simulations we generate 1,000 realiza- 
tions of 5° X 5° convergence maps for each of the three 
source redshifts. It was shown that the ray-tracing simu- 
lations are reliable within a 5% accuracy up to Z ~ 6, 000 
and I ~ 4, 000 at Zg = 1.0 in terms of t he p ower spectrum 
and the bispectrum, respectively [see. [TtI [35j. 

In addition to the fiducial cosmology, we also per- 
form ray-tracing simulations for several cosmologies with 
slightly different parameters. We vary each of the fol- 
lowing cosmological parameters: As, n^, the cold dark 
matter density Q.ch?, f^A, and wq by ±10%, and Wa by 
±0.5. In varying the parameters, we keep the flatness of 
the universe as well as the physical baryon density il^K^ 
unchanged. Therefore, the three parameters, h, r^m, and 
rib, are varied simultaneously to satisfy the above condi- 



tion. For each of these 12 different cosmologies, we gen- 
erate 40 realizations of convergence fields for each of the 
three source redshifts. See Sato et al. [13 for more de- 
tails of the methods used for the ray-tracing simulations 
(see also Sato et al. [26(1). AH the convergence maps used 
in this paper are the same as those used in Seo et al. p^ . 

We include only the auto spectra in our tomo- 
graphic analysis. We assume a future, wide-field weak- 
lensing survey of 1,500 deg^, expected for Subaru Hy- 
per Suprime-Cam (HSC) Weak Lensing Survey [1] for 
the signal-to-noise ratio and the Fisher matrix analyses 
presented below. We simply scale each element of the 
covariance matrix obtained from ray-tracing simulations 
by the ratio of the area. 1/(1,500/25), although strictly 
speaking, the HSV terms have a different scaling. The 
inaccuracy of the above scaling is shown to have little 
impact on the estimation of signal-to-noise ratio (see the 
left panel of Figure 11 in Kayo et al. [iBl)- For cosmologi- 
cal parameter constraints, this difference is much smaller 
than the signal-to- noise ratio as discussed in Section IIXI 
Therefore, we use this simple scaling for the covariance 
matrix, which does not change the results quantitatively. 

In reality, the observed power spectrum is contami- 
nated by the intrinsic ellipticity noise. For the simulated 
convergence map, we can include the noise contamina- 
tion by adding, to each pixel, a random Gaussian noise 
with variance 
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(t) 



(16) 



where is the rms of the intrinsic shear, fig is the 
mean number density of source galaxies, and f2pix is the 
pixel area. Here, we set the intrinsic shape noise as 
Cint = 0.374 and the shear responsivity as 7^ = 1.7 [stI ]. 
We also adopt fig = 12.75, 7.91, and 9.0 arcmin"^ 
at three source redshifts, 2s=0.6, 1.0, and 1.5, respec- 
tively. These values are roughly expected in Subaru Hy- 
per Suprime-Cam Weak Lensing Survey Q and calcu- 
lated from Equation (20) in Takada and Jain [38|, by di- 
viding source galaxies into < Zg < 0.8, 0.8 < Zg < 1.2, 
and 1.2 < z,. 



B. Analysis 

The binned power and bispectrum of the convergence 
field are measured from the simulations as follows. We 
first apply Fast Fourier Transformation to each of the 
convergence fields to obtain k{l). We then bin the data 
into logarithmically-equal bins in /, whose width are set 
as A\nl = In 2/2 « 0.35. The power spectrum and the 
bispectrum of the mth realization are obtained by simply 
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averaging the products of modes: 

^TW^^E 1^(01', (17) 
B':{hMM)^jT^ J2 ^(^{HhMhMh)], (18) 

where Re[...] denotes the real part of a complex number, 
and the summation runs over modes I {U, i = 1,2,3) 
which falls into bin / (li) for the power (bi-) spectrum. 



Cov[PJ0,-Pk(O] = 
Coy[B4lul2,l3),B4l[J'2,l',)] = 



where Nr is number of realizations. Note that these are 
unbiased maximum-likelihood estimators for the covari- 
ance matrices. 

IV. COMPARISON OF FITTING FORMULAE 
WITH RAY-TRACING SIMULATIONS 

We now compare the convergence power and bispec- 
tra measured from our ray-tracing simulation with fitting 
models, which have been used in some previous works. 

A. Weak lensing power spectrum 

Figure [T] shows the convergence power spectrum ob- 
tained from ray-tracing simulations with 1-cr error bars 
expected from a HSC type survey, i.e., fig = 1,500 deg^ 
for Zs ~ 1.0. The numerical error bars increase on large 
scales because of the finite size of the simulation box. We 
compare the simulation result with two fitting formulae 
for the matter power spectrum, revised halof it (solid 
line) and halof it (dashed line). For reference we also 
plot the linear power spectrum result as long-dashed line. 
The vertical arrow indicates the multipole below which 
the simulation result is consistent with higher-resolution 
simulations (512^particles) within 5%. 

The non-linear gravitational evolution of the matter 
density field amplifies the weak lensing power spectrum 
on small scales. We recover a well-known fact that the 
halof it undcrpredicts the convergence power on small 
scales I > 3, 000 [e.g., [H [H [M [H , while the revised 



In the above, Ni is the number of modes taken for the 
summation. Similarly, the factor Ni-^j^^i^, which appears 
in the estimator of the bispectrum, denotes the number 
of triangles in I space. We then average the measured 
spectra over 1, 000 random realizations to obtain our final 
estimates of and P^- We also estimate the the full 
covariance matrix of the power and bispectrum including 
not only the covariance between two different bins inside 
each of the two spectra, but also their cross covariance 
using the 1,000 realizations as follows: 



I 

halof it shows a better match to the simulations. This 
consistency with independent simulations in the litera- 
ture assures the reliability of our simulations. 

Another option to obtain an accurate prediction of 
the power spectrum over a wide range in multipole 
I is to use the model combining perturbation theo- 
ries at small I and halo model at large I, proposed 
by Valageas and Nishimichi [4^ (hereafter combined 
theory). However, it is not straightforward to calculate 
the combined theory, especially its 2-halo term, because 
it involves some renormalization techniques with time- 
consuming multi-dimensional integrals. We thus adopt 
the fitting formulae in this paper for simplicity, but see 
Valageas et al. [Hj for an extensive comparisons with this 
combined theory with the same numerical simulations 
as presented in this paper. 

B. Weak lensing bispectra for equilateral triangles 

Left panel in Figure [2] shows the convergence bispec- 
trum measured from the 1,000 ray-tracing simulations 
at source redshift Zs = 1.0 for equilateral triangles. We 
compare it with several fitting formulae to check the va- 
lidity and usefulness of them. The error bars are ex- 
pected errors for fig = 1,500 deg^ as before. Theoret- 
ical predictions are calculated from Equation ([7]) using 
the fitting formulae by SCOl and Gil-Marinl2 for Bg. 
The solid and dashed curves show the SCOl fitting for- 
mula, while the dot-dashed and dotted curves are the 
Gil-Marinl2 fitting formula with revised halof it and 
original halof it, respectively. We also plot the tree- level 
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FIG. 1: Convergence power spectrum at source redshifts 
Zs = 1.0. The cross symbols are the results of the ray-tracing 
simulations with 1-a error bars assuming the survey area of 
fls = 1,500 deg^. The solid and dashed lines show the two 
fitting formulae, revised halof it and halof it, respectively. 
We also show the linear theory prediction by long-dashed line. 
The vertical arrow denotes the multipole I up to which the 
ray-tracing simulation result is accurate within 5%, which is 
determined based on a convergence test using higher resolu- 
tion simulations. 

perturbation theory (4l| in long-dashed lines. It is shown 
that non-linear gravitational clustering significantly en- 
hances the bispectrum amplitude compared to the tree- 
level perturbation theory prediction by more than an 
order of magnitude at I ^ 500. The vertical arrow at 
/ ^ 4, 000 shows the multipole I up to which the ray- 
tracing simulation result is valid within 5% for the bis- 
pectrum [see, [3^. 

From this figure, it is not obvious to conclude which fit- 
ting formula gives a better prediction, because the power 
spectrum in evaluating the fitting formula of the lens- 
ing bispectrum is also important especially at non-linear 
scales. Therefore, we examine a different quantity which 
we expect to depend on the choice of the power spec- 
trum only weakly in what follows. In order to remove 
the leading, quadratic dependence of the bispectrum on 
the power spectrum, we consider the reduced convergence 
bispectrum defined as 

(-) (, 1 / ^ B^[li,h,h) 

(22) 

which reduces to _B„(Z, Z, Z)/3Pk(/)^ for equilateral trian- 
gles. 

The right panel in Figure [5] shows the reduced con- 
vergence bispectrum Qk measured from the ray-tracing 
simulations averaged over 1,000 realizations at source 



redshift Zs = 1.0 for equilateral triangles. The symbols 
and lines are the same as in the left panel. First, as ex- 
pected, this quantity is much less sensitive to the choice 
of the power spectrum used in the formulae. Our simula- 
tion results are consistent with the Gil-Marinl2 fitting 
formula results up to Z ;$ 4, 000 within the error bars, 
whereas SCOl fitting formula underestimates the ampli- 
tude of the reduced convergence bispectrum for equilat- 
eral triangles, although the broadband shape looks very 
similar to Gil-Marinl2 fitting formula. This is probably 
because Gil-Marinl2 use simulations with more parti- 
cles in larger boxes, and they performed more realiza- 
tions compared with SCOl^, which results in a great im- 
provement in estimation of the bispectrum. Therefore, 
Gil-Marinl2 fitting formula is more suitable to estimate 
the bispectrum and gives larger power than SCOl fitting 
formula. Note also that as shown in Figures 1 and 2 in 
Valageas et al. ^5|, we can also accurately predict the 
reduced bispectrum using the combined theory. 

The power and bispectrum obtained from simulations 
are underestimated at scales I > 6, 000 and I > 4, 000 due 
to the triangular shaped cloud assignment scheme used 
to obtain two-dimensional gravitational potential of the 
lens plane (see, Sato et al. |17|). 

C. Weak lensing bispectra for isosceles triangles 

In Figures [3] and 21 we show the reduced convergence 
bispectrum at source redshift Zs = 1.0 for isosceles 
triangles, where h = l2- As for the equilateral config- 
urations shown in Figure [21 the tree-level perturbation 
theory does not show a good agreement, and we thus do 
not plot it. We also do not plot results of fitting formu- 
lae with halof it, because both the results with revised 
halof it and with halof it are nearly identical. In the 
Figure [31 where ^3 = 152 (left panel), ^3 = 431 (mid- 
dle panel), and Z3 = 1,218 (right panel), we find that 
Gil-Marinl2 fitting formula well reproduces the simu- 
lation results up to li = I2 ^ 1,000, while SCOl fitting 
formula agrees with simulation results up to h = I2 ^ 300 
and underestimates the bispectrum in general. 



^ Gil-Marinl2 used two different simulations named by A and 
B, where a number of realizations are 40 and 3, and combined 
those to obtain a fitting formula. They employ the same ACDM 
cosmology for both set of simulations. Simulations "A" were 
performed with 768^ particles in cubes with 2400/i~^Mpc on a 
side, while "B" adopt 1024'^ particles with the side length of 
1875/i~^Mpc. In contrast, SCOl performed A'^-body simulations 
for various cosmological models, but they performed only one re- 
alization for each cosmological model. Each simulation has 256'^ 
particles in a cubic box 240/i~^Mpc on a side. According to the 
Gil-Marinl2, SCOl fitting formula underestimates simulation re- 
sults up to 20% and biggest discrepancy from the simulations are 
observed at equilateral configuration for both formulae. Consis- 
tently to the previous findings, we can confirm that the predic- 
tion of SCOl underestimates our simulation results as shown in 
Figure [2] 
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FIG. 2: Left and right panel show the convergence bispectrum and reduced bispectrum at source redshift = 1.0 for 
equilateral triangles. The solid and dashed curves are the results of Scoccimarro and Couchman [Sj fitting formula, while 
the dot-dashed and dotted curves are those of Gil-Marm et al. [2^ fitting formula with the revised halofit and original 
halof it, respectively. The tree-level perturbation theory prediction is plotted as long-dashed lines. The vertical arrow shows 
the scale up to which the ray-tracing result is valid within 5%. 
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FIG. 3; The reduced bispectrum as a function of h = I2 at source redshift Zs = 1.0 for isosceles configurations, where Z3 
is chosen to 152 (left panel), 431 (middle panel), and _1,218 (right panel). The solid and dot-dashed lines are the results of 
Scoccimarro and Couchman [S^l and Gil-Marm et al. [23I fitting formula with the revised halofit. 



In the Figure m we show the dependence on for fixed 
h = I2 = 431 (left panel), /i = ^2 = 1, 218 (middle panel), 

= I2 ^ 3,444 (right panel). Again, Gil-Marinl2 fit- 
ting formula is in fairly good agreement with the sim- 
ulation results except for li = I2 = 3,444, while SCOl 
fitting formula underestimates the bispectrum in general 
as in Figure[3l Gil-Marinl2 fitting formula generally de- 



viates from simulation results in squeezed configurations 
h = h ^ 10 ^3. Note that the combined theory provides 
a better match to the simulation results even in squeezed 
configurations It_ = h ~ 10 ^3, as shown in Figure 3 in 
Valageas et al. [35{ . 
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FIG. 4: The reduced bispectrum as a function of /a at source redshift Zs — 1.0 for isosceles configurations, where h — h is 
chosen to 431 (left panel), 1,218 (middle panel), and 3,444 (right panel). The symbols and lines are the same as the Figure|3] 



V. COVARIANCE MATRIX OF WEAK 
LENSING BISPECTRUM 

In this section, we investigate how large the non- 
Gaussian error of the covariance matrix of the lensing 
bispectrum is, compared to the Gaussian error. Wc fo- 
cus on the results of equilateral configurations. See Kayo 
et al. [l6j for the bispectrum covariance matrix of non- 
equilateral triangle configurations. 



A. Diagonal components of the covariance matrix 

Figure [5] shows the diagonal elements of the covariance 
matrix for the convergence bispectrum with and without 
shape noise. We plot the results for equilateral trian- 
gles at Zs = 0.6, 1.0, and 1.5 as a function of multipole. 
The values are divided by the Gaussian contribution of 
the covariance matrix, which is computed by inserting 
the non-linear power spectrum measured from the ray- 
tracing simulations into the first term in the right-hand 
side of Equation ((T^ . Thus, the relative amplitude of the 
non-Gaussian terms is indicated by the deviation from 
unity in this figure. When we neglect the shape noise, wc 
can sec that the non-Gaussian terms become significant 
at multipoles of a few hundreds, and then they dominate 
over the Gaussian component on smaller scales and at 
lower source redshifts due to the non-linear evolution of 
the matter clustering. This trend is similar to that in the 
power spectrum covariance examined in Sato ct al. (iTj 
(see Figure 6 in that paper). 

However, when we add realistic shape noise contami- 
nation described in Section UlI Al the values significantly 
decline and approach to unity irrespective to the redshift. 
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FIG. 5: Diagonal components of the covariance matrix of the 
lensing bispectrum with and without shape noise at source 
redshifts Zs — 0.6, 1.0 and 1.5 for equilateral triangles. The 
results are divided by Gaussian covariances denoted by the 
first term on the right-hand side of Equation H12p . Therefore, 
the deviations from unity arise from the non-Gaussian error 
of the covariance. 



This result is attributed to the fact that the shape noise 
dominates the total error budget, which follows Gaus- 
sian statistics and contributes to the both numerator 
and denominator. Although the importance of the non- 
Gaussian terms arc largely degraded in the presence of 
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the shape noise, we will later show that the impact of the 
non-Gaussian error can remain significant in terms of the 
signal-to-noise ratio as well as the estimated error on the 
cosmological parameters. 

B. Off-diagonal components of the covariance 
matrix 

The correlation coefficient between the convergence 
bispectrum covariances at different triangular shapes 

I 



where li = I'l, I2 = 1'2, and ^3 = Z3 gives the diago- 
nal components, which equal to unity by definition. For 
the off-diagonal components r ~ 1 (— f) means a strong 
(anti-) correlation between the two triangles, while r = 
means the statistical independence of the two. 

Figure [6] shows the correlation coefficient matrices of 
the weak lensing bispectrum for equilateral configura- 
tions at Zs = 0.6, f.O, and f.5. The upper triangular 
parts of the matrices are the results including the shape 
noise, while the lower triangular parts ignore that. In 
the absence of shape noise, the correlation is stronger 
at higher multipolcs and at lower rcdshifts, as expected. 
Compared to the results in the power spectrum case (see 
Figure 8 in Sato et al. the relative strength of off- 

diagonal parts is weaker than that of the power spectrum. 
Strictly speaking, the result depends on the bin widths 
but the above statement should be true, because the bin 
widths are almost the same (AlnZ w 0.35 in this pa- 
per while AlnZ ~ 0.3 in Sato et al. Considering 

I 



2 



and 




I[<l'2<l3<l, 



where Cov^^ is the inverse of the covariance matrix and 
we take account of the bins of the power and bispectrum 
in the range 72 ^ I < ^max {I = 72 is the fundamental 



quantifies the relative strength of the off-diagonal com- 
ponent to the diagonal component. We define the corre- 
lation coefficient as 



(23) 

I 

a realistic shape noise expected in a future weak lensing 
survey, we can clearly see that the non- Gaussian correc- 
tions quickly diminish and the off-diagonal components 
approach to zero. 



VI. SIGNAL-TO-NOISE RATIO 

Although we have shown the impact of the non- 
Gaussian correction to the covariance matrix for each 
clement, it is not clear how important it is to understand 
the whole statistical property of the convergence field. 
One of useful quantities to see this is that the signal-to- 
noise ratio (S/N) of the spectra that quantifies the signif- 
icance of the fluctuation. The S/Ns of tomographic lens- 
ing power spectrum and bispectrum are defined as [e.g., 

a. 



(24) 



(25) 



I 

mode of our ray-tracing simulations. If ~ 2tt/5° ~ 72). 
We impose the condition h < I2 ^ h for the bispec- 
trum so that every triangle configuration is counted just 
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FIG. 6: Correlation coefficient matrices of tlie bispectrum covariance obtained from large number of ray-tracing simulations 
with (upper triangular parts of the matrices) and without shape noise (lower triangular parts) at source redshift Zs — 0.6, 1.0 
and 1.5 for equilateral triangles. 



once. The S/Ns are expected to be independent of the 
bin width, as long as the convergence power and bispec- 
trum do not rapidly vary within bin width, or, in other 
words, the bins arc thin enough. 

We can also define a similar quantity for a joint mea- 
surement of the lensing power and bispectrum tomog- 



raphy. When we consider the non-linear growth of the 
matter density field, it is easy to show that the two spec- 
tra are not independent each other and thus the total 
S/N is not a sum of two S/Ns due to the existence of 
the cross-covariance. We define the data vector for the 
joint measurement as 



, m 

I 



where the indices np and ub arc numbers of bins for the 
power and bispectra. The covariance matrix of the vector 
D can be expressed as 



Cov^ 



CovP" 
CovP'-B. 



CovP''B„ 
CovB" 



(27) 



where Cov^"^" is the cross-covariance between the lens- 
ing power and bispectrum. The S/N for the combined 
measurement is then defined as 



described in Section [ill Bl When the number of indepen- 
dent realizations are not so large compared to the dimen- 
sion of the covariance matrix, the inverse of the resulting 
covariance matrix, which will be used for S/Ns and the 
Fisher matrix analysis, would be overestimated, leading 
to incorrectly tight constraints on the cosmological pa- 
rameters. Therefore we would like to correct this effect 
by multiplying a factor shown in Hartlap et al. [l^l . For 
Nfj independent simulations, an unbiased estimator of 
the inverse covariance is as follows 1421 . 



S_ 

N 



E 



E 



Di 



Gov 



Pk+B„ 
l.V .Zs.z 



Dv 



(28) 



According to Hartlap et al. [4^ , the inverse of an unbi- 
ased estimator of Gov does not yield an unbiased inverse 
covariance matrix, Cov~^, in general, although we have 
adopted an unbiased estimator of the covariance matrix 



Gov I unbiased — 



N. 



1 



(29) 



for Nn — 2>p, where p is the number of bins in the spec- 
tra. For our tomographic analysis with Imax = 2, 000, 
the dimensions of the resulting covariance matrices are 
30 X 30, 345 X 345, and 375 x 375 for the power spec- 
trum, bispectrum, and their joint covariance, respec- 
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FIG. 7: The cumulative signal-to-noise ratios (S/Ns) for the 
convergence power spectrum P^, bispectrum and their 
joint measurement Pk + Bk for tomography with three source 
redshifts. The spectra information over a range of multipoles 
72 < I < Zmax is included. The box, circle and cross symbols 
are the simulation results for the S/Ns of the power spectra, 
the bispectra and the joint measurements of the power spectra 
and bispectra, respectively. The dashed, dotted and solid 
curves show the S/Ns of the power spectra and bispectra and 
the joint measurement for the Gaussian covariance cases. We 
assume f2s ~ 1500 deg'^ for the HSC survey and include the 
shot noise contamination to the covariance matrices. 



tivcly. With our 1,000 independent realizations (Nfj — 
1, 000), the correction factor can be important especially 
when we take the bispectrum into account. Therefore, we 
correct the inverse covariance matrix by using the above 
equation. 

Figure [7] shows the S/Ns for the measurements of the 
power spectra, bispectra, and their joint analysis for a to- 
mographic survey with ftg = 1,500 deg^ area as a func- 
tion of the maximum multipole l^aax- We simply scale 
each component of the covariances to obtain the predic- 
tions for fig = 1,500 deg^, as described in Section UlI Al 
Note that we include the shot noise contamination to the 
covariance matrices in this figure. The box, circle and 
cross symbols are the simulation results for the S/Ns of 
the power spectra, the bispectra and their joint measure- 
ment respectively. The dashed, dotted and solid curves 
show the respective S/Ns when only the Gaussian term is 
included in the covariance matrix. Our simulation results 
show that the true S/Ns with a correct non-Gaussian co- 
variance deviate significantly from that with the Gaus- 
sian covariance. The impact of the non-Gaussian covari- 
ance to the total information content is larger for the 
bispectrum than for the power spectrum when all the tri- 
angle configurations are taken into account. It degrades 



the S/N by a factor of 3 (2) for the bispectrum (power 
spectrum) at I < Zmax = 2, 000. The bispectrum adds 
a new information to the power spectrum and increases 
the value of S/N by about 10% for this Zmax compared 
to the power spectrum measurement alone. Note that 
^max ^ 2, 000 is a typical maximum multipole for the up- 
coming surveys and above that multipole we can not gain 
much information because of the large shot noise contam- 
ination as shown in the figure. The trend that the S/Ns 
do not increase significantly at multipoles I ^ 1000 due 
to the non-Gaussian contribution and shape noise con- 
tamination is similar to the results by Kayo et al. [l6j 
in which the S/N is examined for a single source plane 
without tomography. 



VII. FISHER MATRIX ANALYSIS 

In this section, we employ the Fisher matrix analysis 
to show how much cosmological information we can ex- 
tract from the convergence power and bispectrum with 
a special attention to the impact of non-Gaussian er- 
ror contribution. Although the Fisher matrix forecast 
involves some approximations, Wolz et al. j43j recently 
showed that it would be accurate enough for the predic- 
tions of the cosmological constraints from weak lensing 
surveys because the likelihood function is close to Gaus- 
sian. The Markov Chain Monte Carlo method is more 
robust by direct sampling of the full likelihood function 
without approximations [e.g.. l43 - l46j . In this paper, how- 
ever, we simply propagate the errors on the convergence 
power spectrum and/or bispectrum into projections of 
cosmological parameters using a Fisher information ma- 
trix formalism. 

In the Fisher analysis, we need derivatives of the spec- 
tra with respect to the cosmological parameters as well 
as the covariance matrix. We estimate the derivatives 
by taking the difference of the spectra from ray-tracing 
simulations with different cosmological parameters. We 
first discuss this numerical derivatives before showing the 
results of the Fisher analysis. Figure |S] shows the deriva- 
tives of the convergence power spectrum (left panel) and 
bispectrum (right panel) with respect to six cosmological 
parameters in absolute values at = 1.0. To calcu- 
late the numerical derivatives, we use central difference 
method defined as 



9a:„,2^(p) _ a:k,2^(pq -i- ha) - x^^^Xva - ha) 



dpa 



2K, 



(30) 



where A"„.2^ denotes the convergence power or bispec- 
trum and ha denotes the variations of the cosmological 
parameters pa- In this paper, we choose ha to be 10% 
of the fiducial values of pa, but ha = 0.5 for Wa as de- 
scribed in Section lillAI The solid lines show the predic- 
tions of the fitting formulae (revised halof it for left 
panel and Gil-Marinl2 with revised halof it for right 
panel). The dashed lines are the same as the solid lines, 
but we use the central difference method to calculate the 
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FIG. 8: Derivatives of the convergence power spectrum (left panel) and bispectrum (right panel) with respect to various 
cosmological parameters at Zs = 1.0. The symbols are the results from ray-tracing simulations with different cosmological 
parameters. The solid curves are the results of revised halof it for the power spectrum and Gil-Marinl2 with revised 
halof it for the bispectrum. The dashed curves are the same as the solid curves, but we use the central difference method 
(Equation I30p with the same step sizes ha used for the simulations. 
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FIG. 9: Cosmological parameter constraints from the lensing power spectrum tomography using the non-Gaussian (red) and 
Gaussian (blue) covariance matrices. Left panel: without shape noise. Right panel: with shape noise. We use the power 
spectrum information up to /max = 2, 000 assuming a survey area of l,500deg^. 



derivatives with same ha values used in the simulations in 
order to show the convergence of our numerical deriva- 
tives. Both solid and dashed lines coincide with each 
other, although small differences can be seen in and 
Wa- For both the power and bispectrum, the fitting for- 
mulae show good agreement with simulation results, al- 



though we can see large scatters around the results of 
the fitting formula of bispectrum on large scales proba- 
bly because of the small number of realizations (i.e., 40). 
We will discuss the impact of the step sizes of the central 
difference method and scatters of the bispectrum on cos- 
mological parameter estimation below around Figure 1101 
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We now consider a power spectrum analysis based on 
a tomographic survey with three source redshifts. Using 
the numerical non-linear derivatives and the measured 
covariance matrices, we calculate the Fisher information 
matrix as [e.g., H^l 



Ll'<l 



1 ,Z' 



(31) 

We choose /max = 2, 000 for our analysis. The error 
on the Qfth parameter including marginalization over un- 
certainties in other parameters is estimated as cr(j)a) = 



(F ^)aa; where ^ is the inverse of the Fisher ma- 
trix. 

Figure [9] shows the marginalized error constraints with 
(right panel) and without shape noise (left panel) for 
fig = l,500deg^. The figure shows the marginalized la 
error contours from the weak lensing power spectrum 
tomography for every pairs of six cosmological param- 
eters. The red solid lines show the results obtained when 
we adopt the non-Gaussian covariance matrix measured 
from the simulations while the blue solid lines show those 



with Gaussian covariance matrix. As clearly seen in the 
left panel of Figure El the impact of non- Gaussian errors 
are crucial for constraining the cosmological parameters 
accurately in the absence of shape noise. The presence of 
the non-Gaussian components in the covariance matrix 
enlarge the error ellipses of the dark energy parameters 
(wq and Wa) typically by 40%. However, in the pres- 
ence of shape noise, the two contours arc nearly identical 
and the impact of the non-Gaussian covariance is small. 
The largest difference in the parameter constraint is seen 
in r^A, and is about 10%. It might be enough to con- 
sider the Gaussian component of the covariance matrix 
in estimating the cosmological parameters for ongoing 
surveys, as indicated at the end of Section IVl depending 
on at what accuracy one hopes to constrain the param- 
eters. Although current ongoing surveys seem to be fine 
with the Gaussian approximation in the covariance ma- 
trix, it could be crucial to properly take account of the 
non-Gaussianity for ultimately large survey projects in 
future. 

As in the convergence power spectrum, the Fisher ma- 
trix for the convergence bispectrum tomography is given 
as 



E 



E 



dBf, zs{h,h,h) ^ -1/, , , ,/ ,/ ,x 

— Gov (li,l2,h,Zs,li,l2:h^^s) 



h<l2<h<l„ 
I' <1'<1'<U 



dpa 



dpp 



(32) 



where we have imposed the condition li <l2 < h so that 
every triangle configuration is counted once. 

Figure[Tn]is the same as FigurelHl but the result for the 
weak lensing bispectrum tomography. First, we can see 
that the marginalized error contours for the bispectrum 
analysis are narrower than those for the power spectrum. 
This is quite surprising since the bispectrum is higher 
order spectrum than the power spectrum, so one might 
believe that the constraints from the power spectrum is 
weaker for the bispectrum. Indeed, we have shown that 
the S/N for the bispectrum is smaller than that for the 
power spectrum in the last section (see Figure [71) . The 
smaller error ellipses from the bispectrum analysis may 
be attributed to the fact that the bispectrum is more 
sensitive to the cosmological parameters, since the power 
spectrum and bispectrum are proportional to and 
ri^, in spite of smaller S/N from the bispectrum. A 
consistent result, the bispectrum is more powerful to con- 
strain the cosmological parameters, can be found in Berge 
et al. [13, although they did not assume a tomographic 
survey. We confirmed that this is the case for each of our 
convergence maps at three source redshifts. 

However there arc two concerns: the derivatives of 
the Fisher matrix and the redshift distribution of source 
galaxies. 



• We used only 40 realizations for each of the three 
source redshifts for varied cosmologies to estimate 
the derivatives of the Fisher matrix. As shown 
in Figure [51 the derivatives of the lensing bispec- 
trum estimated in a finite volume have large scatter 
around theoretical predictions on large scales. The 
total volume might not be large enough to converge 
the derivatives with our 40 realizations and this can 
lead to an inaccurate constraints. Furthermore, we 
varied each of the cosmological parameters by ± 
10% except for Wa- From Figure [51 we find that a 
10% change in some of the parameters gives a big 
change in the bispectrum. A 10% change in the 
spectral index n^., especially, results in a ^ 100% 
change in the bispectrum from the fiducial cosmol- 
ogy and this might lead to narrower constraints. 

• We used three (5-function like source redshift dis- 
tribution instead of a realistic continuous redshift 
distribution. To check the impact of this treatment, 
we analytically calculate the power and bispectrum 
for two source distributions. We consider source 
galaxies at exactly Zg = 1.0 and a model redshift 
distribution of galaxies used in Takada and Jain (ssj 
(Equation 20 in their paper) where galaxies have 
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In As Wg n^h^ 



FIG. 10: Same as Figure but result for 

a broad distribution in Zs with the mean rcdshift 
of = l-O- We found that the amphtudes of the 
power and bispectrum in the latter case are smaller 
by 30% (~ 60%) for the power (bi-)spectrum 
and thus the derivatives of the Fisher matrix are 
smaller correspondingly. Note that the elements of 
the covariancc matrices also become smaller when 
a realistic source distribution is considered. There- 
fore, the constraints from the power and bispec- 
trum with a realistic redshift distribution would be 
weaker than the constraints presented here and this 
effect should be larger for the bispectrum. 

Therefore, we guess that the sizes of the marginalized 
contours for the bispectrum are misestimated (underesti- 
mated, probably) in our analysis. However, our interest 
here is in how well the approximation of Gaussianity in 
the covariancc matrix is justified for the bispectrum co- 
variance. The above two possible systcmatics arc more 
related to the derivatives of the bispectrum with respect 
to the cosmological parameters, and their impact on the 
non-Gaussian covariancc matrix would not be significant. 

In the presence of realistic shape noise, the impor- 
tance of the non-Gaussian errors are degraded, similarly 
to what we have shown in the power spectrum analysis. 
Again, one might be able to neglect them for the ongoing 
surveys, while they might be important for future wide 
field surveys. 

We finally calculate the Fisher matrix for a combined 
measurement of the lensing power and bispectrum to- 
mography up to Zmax = 2,000. As in the S/N, the total 
Fisher matrix is not a simple sum of the Fisher matrices 
of the power and bispcctra. By using Equations (j26p and 
(|27p , the Fisher matrix for the combined measurement is 




In As Wg E^oh^ 

the weak lensing bispectrum tomography, 
then defined as 

The cosmological parameter constraints from the joint 
measurement of the lensing power spectrum and bis- 
pectrum with (right panel) and without (left panel) 
shape noise for 0,^ = l,500deg^ are shown in Fig- 
ure [TTJ The red solid lines are the results calculated 
using full non-Gaussian covariancc matrix including the 
cross-covariancc, while the blue solid lines are those using 
the covariancc matrix without the cross-covariance be- 
tween the power and bispcctra. As shown by Kayo et al. 

, the leading contribution to the cross-covariance be- 
tween the two spectra comes from the five point correla- 
tion function, and one might take this as a higher-order 
effect. The impact of the cross-covariance, however, is 
not so small even in the presence of shape noise, and 
the difference is more than 20% for Vlch^ which is well 
constrained from cosmic microwave background experi- 
ments (e.g. WMAP or Planck satellite) though. How- 
ever, for the dark energy parameters such as SIa, wq, 
and Wa the difference is only a few percent. Note that 
the analysis which includes the cross-covariance (i.e., the 
red contours) shows tighter constraints because the cross- 
covariance brings new information. Therefore, when we 
constrain the dark energy parameters from a combined 
measurement of the lensing power and bispcctra, we can 
analyze each spectrum independently and then just add 
information of both spectra to obtain a joint constraint, 
if one needs a precision of 10% in parameter constraints, 
because the impact of the cross-covariancc is only a few 
percent. 
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FIG. 11: Cosmological parameter constraints from a combined measurement of the lensing power spectrum and bispectrum 
tomography using full non-Gaussian covariance (red solid lines) and non-Gaussian covariance without cross-covariance matrix 
(blue solid lines) in the absence (left panel) and presence of the shape noise contamination (right panel), respectively. We use 
the spectrum information up to /max ~ 2,000 and assume survey area is l,500deg^. 



VIII. OPTIMIZING THE DARK ENERGY 
FIGURE-OF-MERIT 

The dark energy Figure-of-Merit (FoM) is often used 
in the hterature to characterize the performance of a sur- 
vey on the dark energy parameter constraints. The dark 
energy FoM is defined as (49| 



FoM 



1 



1 



a{wp)a{wa) ^det(Cov[u;o,Wa]) 



(34) 



where Wp is the equation of state of the dark energy at 
the pivot rcdshift, and Cov[uio, Wa] is the 2x2 submatrix 
of the inverted Fisher matrix, for which we use the same 
method used in the previous section. 

In this section, we examine how FoM scales as a func- 
tion of the mean number density of source galaxies and 
the survey area, under the condition that the total ob- 
servation time is fixed. Following Yamamoto ct al. [soj . 
we express the total survey area fig as 



f Field of View\ ' 



Ttotal 



1-1 X j {tcxv.i "I" ^ohj ) 



where we assume HSC wide survey for 200 days in five 
years with a 1.5 degree field of view and an overhead time 
tohj = 0.3 minutes for each band. We denote by texp.j 
the exposure time for each (the jth) band. We com- 
pute the total observation time as Ttotai = 200 (days) x 
9(hours/days) x 0.7 = 1260 hours, assuming that we take 
an observation for 9 hours per night, and the fine day rate 
of 0.7. We consider the photo-z determination with five 



bands (g, r, i, z, y). We fix the exposure time as 10 (20) 
minutes for g and r bands {z and y bands), but change 



the exposure time for i band <e 



We relate the mean 



number density of source galaxies Ug to the exposure time 
for i band as Hol 



35 



in 



0.44 



30 minutes / 



(36) 



We assume for simplicity that the mean source redshift 
is fixed to 2:m=1.0 regardless of the exposure time. 

Figure [T^] shows the dark energy FoM as a function of 
the mean number density of source galaxies (left panel) 
and the survey area (right panel) , obtained from the con- 
vergence power spectrum P^, bispectrum i?^, and their 
joint measurement + Pk for a tomography with three 
source redshifts under the above condition. We keep the 
fraction of the number density at the three source red- 
shifts the same as in Section fill Al The solid lines show 
the results taking account of the non-Gaussian covari- 
ance while the dashed lines are obtained for the Gaussian 
covariance. The optimal survey design can be different 
from one estimated using a Gaussian covariance, which 
is often done in the literature, when we appropriately 
evaluate the non-Gaussian covariance matrix. As shown 
in the Figure I12[ the value of the mean number density 
(survey area) that gives the maximum FoM is smaller 
(larger) when the non-Gaussian error is properly taken 
into account. For the HSC wide survey, the maximum 
FoM estimated from the joint measurement of the power 
spectrum and bispectrum tomography is given at fig 40 
arcmin"^ (Jls ~ 1200 dcg^) whose values are close to the 
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FIG. 12: The dark energy Figure-of- Merit as a function of the mean number density of source galaxies (left panel) and the 
survey area (right panel), obtained from convergence power spectrum P^, bispectrum B^, and their joint measurement + Bi^ 
for a tomography with three source redshifts assuming Subaru Hyper Suprime-Cam (HSC) wide survey whose total observation 
time is fixed. The solid lines are the results of each spectrum considering non-Gaussian covariance while dashed lines are those 
considering Gaussian covariance. 



planed values (rig ~ 20 arcmin"^, fig ^ 1400 deg^) in 
HSC wide survey (see, HSC white paper and/or SSP 
proposal*). 



IX. CONCLUSION AND DISCUSSION 

We have used the Fisher matrix analysis of a tomo- 
graphic survey with three source redshifts to study how 
well the approximation of Gaussianity is valid for con- 
straining the cosmological parameters, by comparing the 
full covariance matrix obtained from a large number of 
ray-tracing simulations with a Gaussian covariance. Be- 
fore examining the non-Gaussian features of the covari- 
ance matrix, we compared the Icnsing power and bis- 
pectrum obtained from the ray-tracing simulations with 
those from several theoretical fitting formulae. The lens- 
ing power spectrum measured from our ray-tracing sim- 
ulations well reproduces the revised halof it by Taka- 
hashi et al. [2^ based on high-resolution large- volume iV- 
body simulations. For the lensing bispectrum, our ray- 
tracing simulation result shows a good agreement with 



http : //www. slac . Stanford. edu/~ oguri/share/hsc/hsc_whitep: 
http: //www. subarutelescope . org/Science/SACM/Senryaku/HSC. 



the theoretical prediction using a recent fitting formula 
proposed by Gil-Marm et al. [25|, while the fitting for- 
mula proposed by Scoccimarro and Couchman Q gives 
sizable disagreement at multipoles I ^ 1,000. 

We found that the non-Gaussian covariance can be sig- 
nificant at high multipoles in the absence of the shape 
noise contamination, while the shape noise expected from 
a realistic survey adds a large contribution to the Gaus- 
sian part, which degrades the relative impact of the 
non-Gaussian corrections to the covariance. The non- 
Gaussian components of the bispectrum covariance ma- 
trix reduce the cumulative signal-to-noise ratio by a fac- 
tor of 3 at /,jiax = 2, 000 even in the presence of shape 
noise contamination to the covariance. Thus, they are 
crucial to estimate the information content of the lensing 
bispectrum accurately. 

Unlike the S/N, the non-Gaussian terms in the covari- 
ance matrix for the bispectrum degrade the parameter 
constraints by typically less than 15%, which is much 
smaller than that in the S/N. Following Takada and 
Jain [s^, these results can be interpreted as follows. If 
we consider the volume of the Fisher matrix ellipsoid in 
our six-dimensional parameter space as proportional to 
the S/N magnitude, the non-Gaussian covariance shrinks 
the volume by a factor of 3. If all the six principal axes 
,j.(^jjic Fisher matrix ellipsoid are equally shrunk by the 
■opiaait^janssian errors, each parameter error would be de- 
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graded by about 15% 3^^^ — 1), which is close to the 
degradation shown in Figure [TUl When we constrain the 
dark energy parameters such as 51a , wq, and Wa, from a 
joint measurement of the lensing power and bispectrum, 
the impact of the cross-covariance are only a few percent. 

Therefore, the Gaussian approximation of the covari- 
ance matrix might be reasonable for constraining the cos- 
mological parameters in the present weak lensing sur- 
veys. For future weak lensing surveys, however, the non- 
Gaussianity would be considerable when we constrain the 
cosmological parameters with a higher precision. 
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